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Abstract 

The Kardar-Parisi-Zhang (KPZ) equation of nonlinear stochastic growth 
in d dimensions is studied using the mapping onto a system of directed poly- 
mers in a quenched random medium. The polymer problem is renormal- 
ized exactly in a minimally subtracted perturbation expansion about d = 2. 
For the KPZ roughening transition in dimensions d > 2, this renormaliza- 
tion group yields the dynamic exponent z* = 2 and the roughness exponent 
X* = 0, which are exact to all orders in e = {2 — d)/2. The expansion becomes 
singular in d = 4, which is hence identified with the upper critical dimension 
of the KPZ equation. The implications of this perturbation theory for the 
strong-coupling phase are discussed. In particular, it is shown that the cor- 
relation functions and the coupling constant defined in minimal subtraction 
develop an essential singularity at the strong-coupling fixed point. 



1 Introduction 



One focus of today's statistical mechanics is scale invariance far from equilibrium. 
Driven growth of surfaces is an example that widely occurs in nature; for a review, 
see e.g. |]T]]. On large scales of space and time, the effective growth dynamics may 
often be described by a stochastic evolution equation for a continuous "height field" 
h{r,t). The Kardar-Parisi-Zhang (KPZ) equation [Q] 

d^h = uV'^h + ^iVh)'' + r] (1.1) 

driven by Gaussian white noise with 

r)(r,t) = , 

' (1.2) 

r){r,t)ri{r',t') = a^6'^{r - r')6{t - t') 

has become the "standard model" for such processes since it represents the simplest 
universality class of nonlinear growth. Many realistic growth models are in other 
universality classes due to additional symmetries. Moreover, the KPZ equation has 
deep theoretical links with a number of more difficult nonequilibrium problems, 
notably fluid dynamics and turbulence. Phenomenologically, the link to turbulence 
is even more manifest for certain related growth models that show multiscaling . 

The phenomenology of the KPZ equation is well known. In spatial dimensional- 
ities d < 2, the nonlinearity (A/2)(V/i)^ is a relevant perturbation of the Gaussian 
dynamics (A = 0). The strong-coupling regime is characterized by two basic expo- 
nents, the roughness exponent x and the dynamic exponent z, which are defined e.g. 
by the asymptotic scaling on large scales 

{{h{r,, ti) - h{r,, t2)Y) ~ In - r^l'^^Cit \r, - r^in (1-3) 

of the height difference correlation function 0. In the renormalization group, this is 
a crossover between two fixed points: the Gaussian fixed point, which is (infrared-) 
unstable, and the strong- coupling fixed point, which is stable. In ci = 1, the ex- 
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ponents x = and 2; = 3/2 can be obtained exactly in several ways: by ex- 
ploiting the symmetries of the system (namely Galilei invariance and a fluctuation- 
dissipation relation particular to d = 1), by a one-loop dynamic renormalization 
group analysis or by mapping the KPZ dynamics onto an exactly solvable lat- 
tice model P]. All of these tools fail in higher dimensions, and the properties of 
the strong-coupling flxed point are known only numerically. In d = 2, the KPZ 
equation is asymptotically free, and the crossover to the strong-coupling regime is 
exponentially slow 0. Recent numerical values for the exponents are x = 0.386 
and z = 1.612 [0]. For d > 2, the height proflle is smooth in the Gaussian theory. A 
small nonlinearity (A/2)(V/i)^ does not alter this asymptotic scaling; there is now 
a roughening transition to the strong-coupling phase at flnite critical values ±Ac 
[|, ^ |l^. In the renormalization group, the transition is represented by a third 
flxed point. This critical fixed point is unstable and appears between the Gaussian 
flxed point and the strong-coupling flxed point which are now both stable Nu- 
merical studies |Tl|, indicate that a strong-coupling phase with z < 2 persists also 
in high dimensions; various theoretical arguments, on the other hand, predict the 
existence of a flnite upper critical dimension d^, above which 2; = 2 in both the 
weak and strong coupling regimes [|, 0, |l^ . 



A satisfactory theory of stochastic growth should classify the different univer- 
sality classes and the possible crossover phenomena between them, as well as give 
a way to calculate scaling indices exactly or in a controlled approximation. Despite 
considerable efforts, such a theory still seems far. In the framework of the renor- 
malization group, the strong-coupling flxed point does not seem to be accessible 
by the methods of renormalized perturbation theory and the e-expansion that have 
been so successful in equilibrium critical phenomena. This key difficulty is a further 
common feature of the KPZ equation and turbulence, and one may speculate that 
its eventual solution will be similar in both cases as well. 

Recent progress has taken place mainly along a different avenue. The numerical 
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solution of the so-called mode-coupling equations, a self-consistent approximation to 
the full problem, produces exponents which are not very far from the best numerical 



estimates for the KPZ equation |Tj, |T5|, |T^, [T^. Moreover, the mode-coupling 



equations have been shown to be exact in a certain large- limit |15]. However, 
for two reasons it is not clear at present whether the mode-coupling solution can 
serve as the basis for a controlled expansion: (a) The numerical solution of the 
mode-coupling equations presents great difficulties and e.g. for d = 2, it relies so 
far on assumptions on the approximate scaling form of the propagator, (b) The 
transition to finite may well introduce new singularities that have to be handled 
by some kind of renormalization. 

The aim of this paper is to compare the strong-coupling phase with the roughen- 
ing transition under renormalization group aspects. The critical fixed point is per- 
fectly accessible in an e-expansion about the lower critical dimensionality (i< = 2, 
as I show in sections 2 and 3 by dynamic renormalization and by exploiting the 
mapping of the KPZ equation onto a system of directed polymers with quenched 
disorder. The structure of the perturbative singularities is in fact very simple. In 
terms of the dimensionless coupling constant in a minimal subtraction scheme, 
the beta function reads exactly to all orders in perturbation theory 

P{ul) = 2eul, + 2{ul,)' (1.4) 

with e = {2 — d)/2. This yields the dimension-independent critical exponents 

^* = 2, x* = (1.5) 

at the roughening transition, which agree with a scaling argument by Doty and 
Kosterlitz |18[ and with a one-loop dynamic renormalization group calculation 



that has recently been extended to two- loop order ||TP|. 



The perturbation theory for the roughening transition is likely to be exact in 
the interval 2 < d < 4; however, c? = 4 is seen to be a singular point. It is hence 
tempting to identify = 4; see the discussion below. 



In sect. 4, I turn to the consequences of this e-expansion for the crossover to the 
strong-couphng fixed point in d = 2. It proves necessary to carefully distinguish 
between fields and couplings in minimal subtraction and properly renormalized fields 
and couplings defined by their finiteness at a renormalization point. In particular, 
it is shown that the functional dependence of the renormalized coupling constant 
ur and the renormalized height field hpi on their minimal subtraction counterparts, 

ur{um) , hR{hM,UM) , (1-6) 

has an essential singularity at um = 0. (In an ordinary e-expansion, the mapping 
ur{um) is a diffeomorphism of which both the ultraviolet fixed point ur = um = 
and the infrared fixed point u*i^{u\j) are regular points.) This property allows to 
pin down the reasons for the failure of perturbation theory for the strong-coupling 
fixed point. 

The results are summarized and discussed in sect. 5. 



2 The Roughening Transition: Dynamic 
Renormalization 

In this section, I will sketch the dynamic renormalization of the KPZ equation 
[|, ^, |2y, |T9| in a formalism that facilitates comparison with the renormalization for 
the polymer system. 



It is convenient to use the dynamic functional |21 



J VhVhexp - J d'^rdto (^^1 + tho (^/lo - ^V^/iq - y (V/iq)' - p 



{2.1] 



in terms of the field and the "ghost" field h^, which generates response functions, 

/ TV _ N+N \ ^ 5 / \ 

(n*^o(r„to,) n ho{rj,toj)) = l[-—r—^{ J] ho{r,,toj)) . (2.2) 
\i=i j=N+i I i=i \j=N+i I 
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Here the convention has been adopted to absorb all dimensionful constants of the 
linear theory into the "canonical" variables 

to^ut, /lo = (^) ' h, h, (2.3) 

which have dimensions 

o?-2 , d + 2 ,^ 

^0 = 2, -;^o = ^— , + = , (2.4) 

respectively. This convention is standard in field theory, but unfortunately is not 
generally used in the literature on dynamic renormalization, which tends to burden 
the calculations with redundant factors. 

The hnear theory has the response propagator 



Goir2 - n,to2 - hi) = {iho{ri, toi)ho{r2, ^02)) = 

^-(r2-ri)^o 



^(^02 — ^Ol) 

(47r(io2 - ^01))-'^/"° """^^ 



to2 — toi 



(2.5) 



The formal expression for the height-height correlation function 



Co(r-2 - ri, to2 - toi) = {ho{ri, toi)^o(»^2, ^02)) = 

J d'^rdto Go{toi - to, roi - ro)G'o(to2 - ^o, ro2 - vq) (2.6) 

requires for d < 2 (i.e. Xo > 0) the introduction of an infrared cutoff. In a system of 
finite size L with periodic boundary conditions, the stationary correlation function 
at late times ti,t2 is translationally invariant; it has a singularity 

Co(n-r2,toi-to2,L)r.L^^°. (2.7) 

Only the infrared-regularized correlation function remains well-defined in the ther- 
modynamic limit L — > 00: 

Co(^2 - ri, to2 - toi) = lim Co(r2 - n, ^02 - ^oi, L) - Cq{0, 0, L) = 



\r2 - rip^°F(|io2 - toi\/\r2 - . (2.8) 



In the interacting theory, the same subtraction is necessary to define the L-indepen- 
dent stationary two-point function C'{r,to,U()); the higher connected correlation 
functions require infrared regularizations as well. The scale L will also serve to 
generate the renormalization group flow below. 
The canonical coupling constant 

Ao - (^) ^ A (2.9) 

has the dimension e = (o? — 2)/2. We define the dimensionless coupling constant 

uo = XoL' . (2.10) 

The response and correlation functions of the nonlinear theory have the crossover 
scaling form 



\i=i j=N+i 



which can be expressed as the "bare" Callan-Symanzik equation 



:2.ii^ 



dL 



d d d 
+ Y.^j^ + ^0 E^oi^ + /3o(«o)^ ^Xo + N{xo + d) 



N N+N 



X 



UMrj^toj) n ho{rj,toj))^0 (2.12) 

\j=l j=N+l I 



with 

(3q{uq) = LOlUq = euo . (2-13) 

For the infrared-regularized correlators, the explicit dependence on L vanishes in 
the thermodynamic limit. 

In the strong-coupling limit Aq — > oo, these correlation functions develop anoma- 
lous scaling and hence a singular dependence on the bare coupling constant Aq. 
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Renormalization consists in absorbing these singularities into new variables 



hR = Zh{uR)hQ , 

hR = Zj;\uR)ho 

tR = Zt{UR)to , 

Ur = Z{ur)uo . 



[2.U) 



Under this change of variables, Eq. ( |2.12| ) transforms into the renormalized Callan- 
Symanzik equation 



d_ 
'dL 



Ao 



^r,— + z{uR)Y^tRj-— + f5R{uR)- Nx{ur) + N{x{ur) + d) 

• orj . oiRj our 



IN _ N+N 

X ( n ^hR{rj, tRj) n hR{rj, tRj] 

j=N+l 



(2.15) 



\j=i 



with 



d 




euR 




1 - 




ziuR) 


= Zo 


duR 


x{ur) 


= Xo 


duR 



1 



(2.16) 

(2.17) 
(2.18) 



A different but equivalent Callan-Symanzik equation is derived in ref. [113. Notice 
that the renormalization of the ghost field is not independent since the response 
function Gi?(r2 — ri,tR2 — tm,w/j) = {ih{ri,tRi)h{r2,tR2)) always has dimension 



d by its definition. Furthermore, a Ward identity due to Galilei invariance |]19 
enforces the following relation between the Z-factors: 



z = Zi'Z,' 



[2.19) 



By inserting this relation into Eqns. ( [^.16D , ( p^.l7D , and ( [^.1^ ), one obtains 

PiuR) = [-2 + ziuR) + xiuR)]uR (2.20) 
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and hence at any nontrivial fixed point 7^ the exponent identity 

z{u*^) + x{u*r) = 2 . (2.21) 
The renormalized variables (|2.14|) can be defined in a nonperturbative way by 



imposing two independent normahzation conditions e.g. on the infrared-regularized 
correlators in an infinite system, 

L''^'C'j,{mR = {),UR) = L'^'>C',{L,to = 0) , (2.22) 
L''GR{0,tR = L^,UR) = L''Go{0,to = L^) . (2.23) 



L is now an arbitrary normalization scale. An alternative to ( 2.22 ) is the normal- 
ization condition 

L^-'^'idtMiL^UR) = -huR (2.24) 

on the universal finite-size correction to the stationary growth velocity in a system 
of size L, where 6 > is a constant independent of ur that is defined in (|2.26| ) 
below. 

In a perturbative e-expansion, there is an alternative way of constructing the 
Z-factors, namely order by order through a minimal subtraction prescription. This 
makes the rhs. of Eqns. ( p.22|) and ( p.23|) analytic functions of the minimally renor- 



malized coupling constant um with coefficients that remain finite as £ ^ 0. As long 
as the coupling constant is small, this scheme is clearly equivalent to normalization 
conditions. As we shall see in sect. 4, this is no longer the case for large values of 
Ur. I will therefore denote all quantities in the minimal subtraction scheme by the 
subscript M, and reserve the term "renormalized" and the subscript R to quantities 
defined by normalization conditions. 

I will discuss the perturbative renormalization not for the momentum-space 
response and correlation functions in the infinite system as it is customarily done 
but for the position space response function G'o('"5 ^o, ^0) and the stationary finite- 
size amplitude {dtQho){L,uo). The calculation is in close analogy to the polymer 
renormalization group of the next section. 
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The response function has the diagrammatic expansion shown in Fig. 1(a). To 
order Mq, the expansion reads 

^^^^0(0, to = i^', uo) = L^GoiO, to = L') + 0{uy, 4) . (2.25) 

The one-loop diagram does not have a pole at c? = 2 since its short- distance sin- 
gularity cancels with a geometric factor 2 — d ||2D|, As Frey and Tauber have 
shown, constructing the strong-coupling fixed point in d = 1 requires taking into 



account this "hidden" pole of the response function |T9|. However, this finite renor- 
malization can be ignored for the critical fixed point above d = 2, which is the focus 
of this section. 

The expansion for the growth velocity is shown in Fig. 1(b). The tadpole di- 
agram at order uq consists of a nonuniversal ultraviolet-divergent part and of the 
universal finite-size correction 

/ d^r'dt', [dr>Go{r'X,L)]' = -b L-" , (2.26) 

J reg. 

At order Mq, the boxed subdiagram contributes a pole, 

L^-^''{dt,ho)iL,no) = -buo (l + -^ui^ + 0{ule'',ul) (2.27) 

with c = l/327r. This pole originates from the integration region where the two 
vertices approach each other. It can be absorbed into the definition of the variables 
hpf = ZMhho, hM = Zj^lyJiQ, tu = ZMtto, um = ZmUq, with 

ZuhiuM) = l-^<, + 0«), (2.28) 
ZMtiuM) = l + 0«), (2.29) 
ZmM = l + ^< + 0«,). (2.30) 



This reparametrization respects ( p.l9|) and renders both the response function and 
the growth rate regular as £ — >• 0, 

L'^Gm(0, tM = L\ um) = L'^GoiO, to = L^) + 0«5°, <) , (2.31) 



L'-^^{dt^,hM){L,UM) = -huM + 0«£°,O 
iFiom. ( p.28|) and (|2.16| ), one obtains the beta function 



Pm{um) = euM + cull + 0{u\i) 



For d > 2, it has a pair of real- valued unstable fixed points 



^2.32) 



(2.33) 



u 



* 2 

M 



(2.34) 



that describe the roughening transition from the weak coupling to the strong cou- 
pling phase. Relations ( p.l7|) and ( p.l8| ) then give the critical exponents 



= 2 + Oie' 



X* = + 0(6' 



(2.35) 



satisfying ( 2.21| ). While higher-order calculations are cumbersome in the dynamic 
framework, we will see in the next section that these values are exact to all orders 
in perturbation theory. 



3 The Roughening Transition: Replica 
Renormalization 

It is well known that the KPZ equation can be mapped onto a system of directed 
polymers given by the partition function 

2 



T>r exp 



(3.1) 



Here r{t) denotes the polymer displacement field in d transversal dimensions as a 
function of the longitudinal "timelike" coordinate t. The polymer is subject to the 
quenched random potential Xri{r,t) that has the statistics ( |1.2|) and describes point 
impurities with short-ranged correlations. 
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The height field of Eq. ( |1.1| ) is related to the restricted partition sum Z{r,t) of 
all paths ending at a given point (r, t) by the Hopf-Cole transformation 



h{r,t) = y logZ(r,t) . 
Hence the disorder- averaged free energy per unit longitudinal length 



(3.2) 



f{L) = -2v lim drlog Z{T,L) 

T — >oo 



(3.3) 



in a system of size TxL (with periodic boundary conditions in transversal direction) 
is proportional to the stationary growth velocity, 



(3.4) 



A convenient way to set up perturbation theory is the replicated partition func- 



tion 



(3.5) 



where the the parameters v and cr^ are again absorbed into the definition of the 
canonical variables to and h^. Since we are interested in arbitrary particle numbers 
p, we rewrite the partition function in second quantization, 



I)(f)V(j) exp 



J dtod'r (0(9,,- 9,2)0 -A^ 



2 12 i2 



(3.6) 



where 0(r, tg) is a complex field. The normal-ordered interaction term —A 



2j,2j,2 ; 



IS 



an attractive pair contact potential. More generally, we define the normal-ordered 
m-line contact fields 



$„(to)= /d'^r[0(r,t)]-[0(r,t)r . 



(3.7) 



The perturbation series for the free energy 

oo „,2N 



U{L, ul) = ^L-'"^' Jdt2... dt^{<l>,m,{t,) . . . $2(t^)), (3.8) 
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is a sum involving connected pair field correlations in the p-line sector of the unper- 
turbed theory {uq = 0). The integrals in Eq. ( ^78[) are infrared-regularized by the 
system size L; their ultraviolet singularities are determined by the short-distance 
structure of the pair field correlations and have to be absorbed into the coupling 
constant renormalization. Hence consider the asymptotic scaling of the A^-point 
function ($2(^01) • • • ^2(toN)) as the points toi, . . . , toAf approach each other, 



tTjk and (3.9) 



with Tjk and the "center of mass" t'^ = N ^ ^f=i ^oj remaining fixed. This is given 
by the p-independent short-distance algebra 

N+l 

$2(toi) . . . $2(W) = E [C™(ri, . . . , r^_2) ^UQ + •••], (3.10) 

m=2 

where are scaling functions of the — 2 linearly independent distance ratios Tjk 
and the dots denote terms that are subleading by positive integer powers of t/L'^. 
The integration over the relative distances then yields 

N-2 

dtt^-^ n dT-i ($2(t0l)...$2(M) = 

1=1 
N+l 

E / 3+-(^-™+i)dt {^miQ) + ■■■ (3.11) 

. n J 



m=2 ■ 



with 



„ N-2 

Jn= li dnC^{n,...,TM-2). (3.12) 

1=1 



Hence we obtain to one-loop order 

L'f,{L,ul) = -L'{<^,), .^ul(l + ^ul)+0iuy,4) (3.13) 



with L'^($2)p = p{p — l)/2. The pole in ( p.l3|) originates from the term in ( |3.11 



with m = N = 2, which is shown diagrammatically in Fig. 2(a). This pole can be 
absorbed into the definition u\.j = Zf.jU^ with Zm given by ( p.30|) . To this order, we 
hence recover the beta function ( |2.33| ) and the fixed points ( |2.34| ) of the dynamic 
calculation. 
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In the polymer framework, however, it is not difficult to discuss higher orders. 
The ultraviolet singularities of the A^-th order integral ( p.ll|) are contained in the 



coefficients or arise from the integration over t. In the ffist case, they are due 
to a proper suhdiagram and hence already absorbed into the renormalized coupling 
constant at lower order. Only the divergences from the integration over t, with J"^ 
denoting the regular part of (|3.12| ), may contribute to the primitive singularity at 
order A^. Inspection of then shows that a pole at e = only appears for 



m = 2. However, there is only one diagram per order of this kind, which is shown 
in Fig. 2(b). Since this diagram factorizes into loops of the kind of Fig. 2(a), it 
contributes a pole in e of order — 1. Therefore the pole at order = 2 is the 
only primitive singularity in the series (|3.8|) for the free energy; analogous arguments 



apply to the expansions of the contact field correlation functions. It follows that the 
one-loop equations ( p. 301) , ( p.33|) and (|2.34| ) are exact to all orders in perturbation 



theory (see also |23|, |24 

The replica trick is unproblematic within perturbation theory, since it reduces 
to convenient bookkeeping of the averaging over disorder. Indeed, the random limit 
p ^ is trivial in Eqns. ( |2.3(J| ), ( |2.33| ) and ( |2.34 ) which are independent of p. The 



crossover scaling function of the disorder-averaged free energy 

C«) ^ L7(L, ul,) = hm -/p(L, O (3.14) 

p^O p 

is a regular function of the minimally subtracted coupling constant u\,j, as follows 
by inserting (g^^) into (|3l3|) . By (jSj), ( ^1281) , (ICTj ), and (g30|), the function 



-C{u\[) /um equals the scaling function ( |2.32|) of the minimally subtracted growth 



rate, and hence 6 = 1/8. If the renormalization point condition ( p.24| ) is chosen, the 
universal function C{u\j) is directly related to the renormalized coupling constant, 
^(■^m) = ^'"/{('"m)- Its finite fixed point value 

= C{u\j) = --^e + 0{e^) , (3.15) 

the Casimir amplitude at the roughening transition, is the analogon of the central 
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charge in conformally invariant field theories [^. The function C{u1,j) has the finite 
hmit (|3.15|) since the free energy does not develop an anomalous dimension at the 
roughening transition, i.e. hyperscaling is preserved. By ( p.4|) , this implies that 
(dth) ~ L^^ at the transition, and hence ~ ^* = 2. The exponents (1^) then 
follow from this relation together with ( p.21|) . 



The dynamic exponent z* = 2 can also be verified independently. Consider the 
two-point function of the normal-ordered density field $i(to) = 00(^ = 0, to), 

oo 2N . 

Y: J^h-.L-'"'^ / dtoi . . . dW($i(to)$i(t[,)$2(toi) . . . $2(W)) . (3.16) 

N=0 ^ ■ 

Its short-distance asymptotics gives the return probability of a single line to the 
origin r = 0. In the linear theory, 

($i(to)$i(t'o))(^,0)~|to-tor'/'° (3.17) 

for |to — ^ol/-^^ ^ 1. Any anomalous contribution to this exponent arises from the 
renormalization of the fields $i(to) and $i(tQ). The renormalization of $i(to) is 
due to a short-distance couphng of the form 

$i(to)<^>2(toi) • • • MtoN) = t-^'^/'CUn, r^_i) $i(to) + . . . (3.18) 

for t/L^ (with toj -to = troj , toj -tok = trjk for j,k = 1, ... ,N, and n, . . . , tn-i 
denoting a basis of the fixed ratios tqj, Tjk), and there is a corresponding expression 
for $i(to). However, it is obvious that the product on the l.h.s. couples only to 
contact fields of at least two lines, and therefore Cj^ = at all orders A^. 



4 The Strong-coupling Fixed Point 

As discussed in the previous sections, the crossover from the critical fixed point 
to the Gaussian fixed point in d > 2 can be parametrized in terms of the coupling 
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constant mm of a minimal subtraction scheme. In the framework of the e-expansion, 
um is completely equivalent to the coupling constant ur defined by normalization 
conditions, e.g. (p.22|) and ( p.23|) : the two couplings are related by a diffeomorphism 



that remains regular in the limit e — > 0, defined on a domain of interaction space 
that contains the fixed points um = ur = and u*^{u\.j). This equivalence is lost 
in the crossover to the strong-coupling fixed point, as I will now show. 

Integration of the flow equation ( |1.4|) for in two dimensions with the initial 
condition u\j{Lq) = ul yields 



l-uf log(L/Lo) 

an expression that diverges at a finite value L = l/uf. It follows immediately that 
um{L) is not well suited to describe the crossover to the strong-coupling fixed point 
since it divides the crossover into an ultraviolet regime L < l/u\ and an infrared 
regime L > l/u\ where u\j{L) is defined by analytic continuation of the solution 
( [4.1| ). However, the pole of u\.j{L) a.t L = l/ufis only a "coordinate singularity" |^ 
of the minimal subtraction scheme; any renormalized coupling ur{L) is expected 
to remain regular at L = Hence the function also has a pole at the 

value u\{L = l/u\) between the fixed points m|. = and u*^^. 

For small u\i > 0, the correlation functions (|2.22| ) and (|2.23|) at the normaliza- 



tion point can be calculated perturbatively, which results in an analytic renormal- 
ization 

hR = [1 + aiu% + a2u'^R + . . .]hM , (4.2) 
tR = [l+b^ul + b2u\+ ...]tM , (4.3) 

and hence by ( |2.21| ) 

uli = [1 + (ai + bi)uR + (^2 + CLibi + b2)u\ + . . . (4.4) 
with finite coefficients a^, b^. Hence has the beta function 

PrK) = 4 + (a? + + bl~a2- b2)ul + 0«) . (4.5) 

15 



In perturbation theory, one would calculate these power series up to some finite 
order in m|, and look for a fixed point of Eq. ( [4.5|) . However, a low-order 
calculation does not yield nontrivial exponents at the strong-coupling fixed point. 
Consider e.g. the expression ( |2.18| ) for the roughness exponent in = 2, 



(4.6) 



It is zero if Zh{u'j^) is a regular function of since (3r{u*^) = 0. In an ordinary 
e-expansion, the beta function and the Z-factor are treated as power series; a finite 
result would then arise from the hnear part of /5r(-u|.) together with the simple 
poles of Zh{uj^). Here, if we assume the existence of an expansion parameter e and 
naively take m|j ~ e, we obtain x{u*ii) = 0(£^°) from (^), (|45|), and (|4^ ). In 



fact, we expect a perturbative fixed point of Eq. (|4.5|) to be spurious at any order 



i.e. not to refiect properties of the large-distance asymptotic regime. The reason 
is that all coefficients ajy, bjy are finite for e = and thus depend on details of 
the infrared regularization, unlike the residues of the poles in e which determine 
the Z-factors in an ordinary e-expansion. It is difficult to see how this dependence 
could cancel out to produce universal exponents. 

In the infrared regime L > Lq, the relationship between u\j and ufj is no longer 
given by Eq. ( [4.4|) . As L —>■ oo, the asymptotic behavior of the renormalized 
quantities is 

<2-4~L-^'=exp[i/7<], (4.7) 

tn = ZttM ~ = exp[-(z - ^o)/mm] , (4-8) 

hR = {uijU'^Z^;^)hM ~ (n^/exp[(2; - ZQ)/ul,^]hM , (4.9) 

in terms of their minimal subtraction counterparts (the exponent y' is defined by 
Pr^'^r) = + Hence all of these quantities have essential 

singularities at um = 0, which are tied to an essential singularity 

{hM{ri,t)hMir2,t)) ~ uj/exp[-2(z - Zo)/ulj] (4.10) 
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of the correlation function in minimal subtraction in the same limit, for fixed dis- 
tance |ri — r2|. This shows the inequivalence of the two renormalization schemes. 

5 Discussion 

In this paper, it has been shown that the perturbation theory for the KPZ rough- 
ening transition can be developed consistently in the dynamic and in the directed 
polymer framework. To all orders, it predicts the exponents 

z* = 2, C = (5.1) 

resulting from the exact beta function (|1.4|) in minimal subtraction. 

This perturbation theory gives some indications on the existence of an upper 
critical dimension for the KPZ equation. At the critical fixed point u\f^ the 
nonlinearity is a relevant perturbation 

= -2e« - O + . . . , (5.2) 

and the interface roughness changes to leading order in perturbation theory, 

x(«m) = « - O + • • • • (5.3) 

With the assumption that x{L) is a monotonic function of L, this implies x > 
in the strong-coupling phase. This argument, however, can only be trusted for 
d < A. The minimally subtracted pair contact field $2,m has the exact dimension 
2 + 2e at the critical fixed point. This dimension would turn negative for d > 4, 
which is obviously impossible. The regularization of the perturbation series ( |3.8| ) 
by minimal subtraction breaks down at = 4 since already the two-loop integral 
/ dto2('^'2(^oi)'^2(^o2)) develops a new pole that is not absorbed into the coupling 
constant renormalization ( p.3(]| ). This exhibits the singular role of dimension 4 and 
suggests the identification d^ = 4, in agreement with recent mode-coupling results 
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1^ . It should be possible to corroborate these results by extending the analysis of 



the roughening transition to dimensions d > 4. 

The form (|1.4|) of the flow equation has an interesting consequence for the desrip- 
tion of the crossover to strong-coupling behavior in ci = 2. The minimal subtraction 
coupling constant u\j{u']^ develops a pole at a finite value u\ = u\ < u*^ of the 
renormalized coupling constant, a singularity not taken into account by existing 
renormalization group treatments of the KPZ equation. That singularity divides 
the crossover into an ultraviolet regime u\ < u\ and an infrared regime u\ > u\. 
Only the former is accessible in perturbation theory: even the summation of an 
infinite number of terms e.g. in the series ( |4.2| ) and ( |4.3| ) can only lead to a relation 
u\{u\j) for u\[ > 0. Whether information on the infrared regime can be gained by 
inverting this relation and then analytically continuing it to values uj^ > u\ remains 
to be seen. 

The precise form ( p..4|) of the beta function is clearly due to the fact that the 
intrinsic dimension of the polymer, or the dimension of time, is one. The disorder 
problem, however, has a natural generalization from directed lines to D-dimensional 
directed surfaces. For D ^ 1, higher-order perturbative singularities are expected, 
which generate further terms in (|1.4|) . This generalization has an interest in its 
own right. It is not clear, however, if it is useful to solve the original problem: any 
perturbative fixed point will tend to infinity as D ^ 1, but the exponents may still 
be well-behaved in that limit. 

Useful discussions with H. Kinzelbach are gratefully acknowledged. 



18 



References 

[1] J. Krug and H. Spohn, in Solids Far From Equilibrium: Growth, Morphology 
and Defects, ed. C. Godreche (Cambridge, 1990). 

[2] M. Kardar, G. Parisi, and Y.-C. Zhang, Phys. Rev. Lett. 56 (1986), 889. 

[3] J. Krug, Phys. Rev. Lett. 72 (1994), 2907. 

[4] D. Forster, D.R. Nelson, M. Stephen, Phys. Rev. A 16 (1977), 732. 

[5] L.-H. Gwa and H. Spohn, Phys. Rev. Lett. 68 (1992), 725. 

[6] T. Nattermann and L.-H. Tang, Phys. Rev. A 45 (1992), 7156. 

[7] T. Ala-Nissila et al., J. Stat. Phys. 72 (1993), 207. 

[8] J.Z. Imbrie and T. Spencer, J. Stat. 52 (1988), 609. 

[9] J. Cook and B. Derrida, Europhys. Lett. 10 (1989), 195; J. Phys. A 23 (1990), 
1523. 

[10] M.R. Evans and B. Derrida, J. Stat. Phys. 69 (1992), 427. 

[11] J.M. Kim and J.M. Kosterhtz, Phys. Rev. Lett. 62 (1989), 2289; B.M. Forrest 
and L.-H. Tang, Phys. Rev. Lett. 64 (1990), 1405; J.M. Kim, J.M. Kosterhtz, 
and T. Ala-Nissila, J. Phys. A 24 (1991), 5569; J.M. Kim, M.A. Moore, and 
A.J. Bray, Phys. Rev. A 44 (1991), 2345; L.-H. Tang, B.M. Forrest, and D.E. 
Wolf, Phys. Rev. A 45 (1992), 7162. 

[12] T. Halpin-Healey, Phys. Rev. A 42 (1990), 711. 

[13] M Feigelman et al., Phys. Rev. Lett. 63 (1989), 2303. 

[14] A. Bouchaud and M. Gates, Phys. Rev. E 47 (1993), R1455; Phys. Rev. E 48 
(1993), 635. 

19 



[15] J.P. Doherty et al, Phys. Rev. Lett. 72 (1994), 2041. 

[16] Y. Tu, Phys. Rev. Lett. 73 (1994), 3109. 

[17] M.A. Moore, T. Blum, and J.P. Doherty preprint (1994). 

[18] C.A. Doty and J.M. Kosterhtz, Phys. Rev. Lett. 69 (1992), 1979. 

[19] E. Prey and U.C. Tauber, Phys. Rev. E 50 (1994), 1024. 

[20] T. Sun and M. Phschke, Phys. Rev. E 49 (1994), 5046. 

[21] P.C. Martin, E.D. Siggia, H.A. Rose, Phys. Rev. A8 (1973) 423; R. Bausch, 
H.K. Janssen, H. Wagner, Z. Phys. B24 (1976), 113. 

[22] H.W.J. B15tc, J.L. Cardy, and M.P. Nightingale, Phys. Rev. Lett. 56 (1986), 
742; I. Affleck, Phys. Rev. Lett. 56 (1986), 746. 

[23] S.M. Bhattacharjee and J.J. Rajasekaran, Phys. Rev. A 44 (1991), 6202. 

[24] M. Lassig, Phys. Rev. Lett. 73 (1994), 561. 

[25] M. Lassig, Nucl. Phys. B 334 (1990), 652. 



20 



Figure Captions 

1. Diagrammatic expansions generated by the dynamic functional ( p.l|) . The 
lines with one and two arrows denote the unperturbed response function (|2.5| ) 
and the unperturbed correlation function ( p.8[ ), respectively. Dots represent 
the vertices ih{WhY\ each incoming line to a vertex has to be differentiated 
with respect to r. (a) Response function Go (r, to, Mo)- The one-loop diagram 
is regular at e = 0. (b) Stationary growth rate {dtohQ){L,uo). The boxed 
sub diagram contains a simple pole at £ = 0. 

2. Singular diagrams contributing to the finite-size free energy fp{L, Uq) in the ex- 
pansion ( |3.8| ). The lines denote unperturbed single-line propagators 
{(f>{ri,toi)(f){r2,tQ2)), the dots pair contact vertices ^2(to)- (a) One-loop di- 
agram containing a primitive pole at e = 0. (b) A^-loop diagram containing a 
pole of order — 1 at e = 0. 
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